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ABSTRACT 



We empirically test the relation between the SFR(L/^) derived from the infrared lumi- 
nosity, Lir, and the SFR(Ha) derived from the Ha emission line luminosity using simple 
conversion relations. We use a sample of 474 galaxies at z = 0.06 — 0.46 with both Ha de- 
tection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 /xm). We 



* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA. 
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derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very 
clear trend between E(B — V) and Lju that allows to estimate extinction values for each 
galaxy even if the H/3 emission line measurement is not reliable. We calculate the L/ r by in- 
tegrating from 8 up to 1000 /im the SED that is best fitting our data. We compare the SFR(Ha) 
with the SFR(Lja). We find a very good agreement between the two SFR estimates, with a 
slope of to = 1.01 ± 0.03 in the log SFR(L/#) vs log SFR(Ha) diagram, a normalization 
constant of a = —0.08 ± 0.03 and a dispersion of a — 0.28 dex. We study the effect of some 
intrinsic properties of the galaxies in the SFR(L/#)-SFR(Ha) relation, such as the redshift, the 
mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR 
comparison. The mean ratio of the two SFR estimators log[SFR(Ljjj)/SFR(HaO] varies by 
~ 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(0/H) + 12 < 9.2). This effect 
is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. 
Considering different morphological types, we find a very good agreement between the two 
SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and 
normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much 
steeper (to = 1.62 ± 0.43), but it is still consistent with 1 at the 1.5a level, because of the 
reduced statistics of this sub-sample. 

Key words: galaxies: formation — infrared: galaxies — (ISM): dust, extinction 



1 INTRODUCTION 

The star formation rate (SFR) is one of the most important parame- 
ters to study galaxy evolution, as it accounts for the number of stars 
being formed in a galaxy per unit time. Many authors have studied 
the evolution of the star formation rate density with time, showing 
that the peak of the star formation took place at z ~ 1 — 2, then 
slowing do wn at later times (see the d ata collection from different 
surveys by iHopkins & Beacomll2006l and recent Herschel results 



surveys by Hopkins & tseacom 2UUb ana recent Herschel results 
from iGruppioni et al.l 201ol) . There is also agreement in the fact 
that SFR directly correlates with the galaxy properties such as gas 



conten t, mo rphology or m a ss. Recent papers s u ch as Elbaz et al 



201 lh . IWuvtsetalJfeOllh . iDaddi et alj hold) . iRodighiero etal . 
20101) or iNoeskd j2009n " have investigated in detail the existence 



of a main sequence in the SFR-mass relation. Therefore, it is fun- 
damental to have a rigorous way to measure SFR when studying 
galaxy evolution. 

A great effort has been made over the past years to derive SFR 
indicators from luminosities at different wavelengths, spanning 
from the UV, where the recently formed massive stars emit the bulk 
of their energy; to the infrared, where the dust-reprocessed light 
from those stars emerges; to the radio, which traces supernova ac- 



Kennicutl 1998, Yun et al. 2001, Kewlev et al. 2002, Ranalli et al. 
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2006, Alonso-Herrero et al. 2006, Rosa-Gonzalez et al. 2007, 
Calzetti & Kennicutt 2009, Cal zetti et alj l20ld) . Recent works 



combine optical and infrared observations to derive attenuation- 
corrected Ha and UV continuum luminosities of galaxies by 
combining these fluxes wit h various components of IR trac 
ers (e.g. 



Bell et al 



Gordon et al.l l2000l llnouel l200lL iHirashita et al 



2003, 



Hirashita et al 
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The SFR derived from the Ha emission line is calibrated on 
the physical basis of photoionization. The nebular lines effectively 
re-emit the integrated stellar luminosity of galaxies short-ward the 
Lyman limit, so they provide a direct, sensitive probe of the young 
massive population. Only stars with masses > lOM© and lifetimes 
< 20 Myr contribute significantly to the integrated ionizing flux, 
so the emission lines provide a nearly instantaneous measure of the 



SFR, independent of the previous star formation history. However, 
a large part of the light emitted by young stars, which mostly re- 
side within or behind clouds of gas and dust, is absorbed by dust 
and then re-emitted at longer wavelengths. To properly calculate 
the SFR from Ha emission line an extinction correction must be 
applied to take this effect into account. 

The absorbed light from these young stars, tracers of the SFR, 
is re-emitted at IR wavelengths. Therefore, the SFR can also be 
derived from the IR emission of the galaxies. One of the most 
commonly used met hods to derive th e SFR from the IR emis- 
sion is to apply the iKennicuttl ( 1998) [K98 hereafter] relation, 
which links the SFR to the IR luminosity (Lie). However, this 
is a theoretical relation base d on the starburst synthesis models of 
iLeitherer & Heckma nldl995h obtained for a continuous burst, so- 
lar abundances and a Salpete n Jl955l) Initial Mass Function (IMF), 
assuming that young stars dominate the radiation field throughout 
the UV- visible and that the FIR luminosity measures the bolometric 
luminosity of the starburst. This physical situation holds in dense 
circumnuclear starbursts that power many IR-luminous galaxies. 
In the disk of normal galaxies or early type galaxies, the situation 
is much more complex as dust heating from the radiation field of 
older stars may be very important. Strictly speaking, as remarked 
by K98, the relation above applies only to starburst with ages less 
than 10 8 years while, for other cases, it would be probably better to 
rely on an empirical calibration of SFR/L/.R. 

The advent of the Herschel Space Telescope, which samples 
the critical far-infrared peak of nearby galaxies, is an excellent op- 
portunity to empirically derive the relation between the observed 
L,iH and the SFR without associated uncertainties due to extrap- 
olation of the IR peak. In this paper we analyze a far-IR (100 or 
160 /im) selected sample in the COSMOS field with associated 
multi-band photometry and 20k zCOSMOS optical spectroscopy 
at z = 0.1 — 0.46. We determine the SFR estimates from the Ha 
line luminosity corrected for extinction, while we obtain the Lj_r by 
integrating from 8 to 1000 /im the best fitting SED of our galaxies. 
We test if the SFR derived from the Ha emission and the SFR de- 
rived from the Lir are consistent at different redshifts and physical 
properties of the galaxies (mass, metallicity, specific star formation 
rate or morphological type). 

This paper is organized as follows: Section 2 introduces our 
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sample of galaxies. Sections 3 and 4 explain the method used to 
derive the fundamental physical parameters of our sample (Lir, 
stellar masses, SFR from dust extinction corrected Ha emission 
line). Section 5 presents the comparison of the two SFR indicators, 
while in Section 6 we discuss the effect of different galaxy proper- 
ties on the comparison of the two SFR estimates. Finally, in Section 
7 we summarize our conclusions. 

Throughout this paper we use a standard cosmology (fl m = 
0.3, Q,a = 0.7), with H = 70 km s _1 Mpc" 1 . The stellar masses 
are given in units of solar masses (M©), and both the SFRs and the 
stellar masses assume a Salpeter IMF. 



2 DATA 

We study the SFR derived from the Lja making use of the re- 
cent Herschel data in the COSMOS field. The Herschel data in- 
cludes PACS jPoglitsch etai] l2010h data fro m the PEP survey 
(PACS Evolutionary Probe. iLutz et al.ll201ll) . We use the PACS 
blind catalog v2.1 (containing 7313 sources) down to 3a, cor- 
responding to ~ 5 and 10.2 mjy at 100 and 160 /im, respec- 
tively. This catalog was associated u sing the likelihood ratio tech- 
nique ^Sutherland & Saunders! [l992h to the deep 24 //m catalog 
from Le Floch et al. (private communicat ion Sii m (24 /im)=80 /iJy; 
iLe Floc'h et al.l2009USanders et al.l2007l) and to the multi-band 3.6 
^m selected catalog from lllbert et al.l d2010h (Szi m (3.6 /_m)=1.2 
fjjy), including NUV, u* , Bj, g+, V Jt r+, i+, z+, J,K S , 3.6, 
4.5, 5.8, 8.0 and 24 from GALEX, Subaru, IRAC and Spitzer 
legacy surveys, respectively. We also need the spectroscopic in- 
formation to accurately measure the SFR from the Ha emission 
line. We have 1654 sources (~ 23 % of PEP catalog) with spectro- 
scopic informatio n from the zCOSMOS 20k spectroscopic survey 
dLillv et alj[2009h . 717 of them with z < 0.46. This cut in redshift 
is necessary to properly observe the Ha emission line (6562.8 A). 
We did not consider sources with no accurate redshifts (flag z ^ 
2.1. lLillvet al]|2009l . ~ 5%), neither sources with S/N < 3 in the 
Ha flux (~ 25%). We did not consider galaxies with very high val- 
ues of aperture correction (> 16; 4 galaxies) to exclude significant 
uncertainties associated to this correction. 

To avoid AGN contamination we have cleaned our sample of 
possible AGN candidates. Type-1 AGN's are easily recognizable 
by their broad permitted emission lines and have been excluded 
originally from the zCOSMOS galaxy sample. We have eliminated 
sources classified as Seyfert2 or LINERS using the BPT diagnostic 



Bongiorno et all d2010l) . as well as XMM-Newton 



120071) and Chandra detected sources dElvis et al.l 
2009). We have also performed a SE P fitting from t he NU V to 
the PEP bands to our sources using the iPolletta et al.l ( 2007) tem- 
plates which include different SEDs, from elliptical and spiral nor- 
mal galaxies, to starburst, composi te, Seyfertl and Seyfe rt2. We 
also used the modified templates by Gruppi oni et al.l d2010h . which 
include a higher IR bump to better fit the observed far-IR data. We 
have carefully inspected the resulting SED fits and we have elim- 
inated those sources with a typical SED of a Seyfert. This proce- 
dure allows us to identify AGN sources which could be missed by 
the other methods (e.g., galaxies with too low S/N to be used in 
the BPT diagrams). We identify 47 AGN candidates (~ 6%) and 
eliminate them. The final sample consists of 474 PEP detected Ha- 
emitting galaxies with 0.06 < z < 0.46. In Table[T]we summa- 



rize some of the main statistical properties of our sample (z, stellar 
mass, SFR, SSFR, metallicty and Ha, H/3 equivalent widthsfl 



3 SED FITTING: INFRARED LUMINOSITIES AND 

MASSES 

3.1 Infrared luminosities 

To derive the Lir value of each s ource we have performed a SED- 
fitting using the Le Phare code dArnouts et al.l l200lL lllbert et all 
2006]), which separately fits the stellar part of the spectrum with 
a stellar library and the IR part (from 7 /im) with IR libraries. 
To perform the SED-fitting we used 4 different infrared libraries 
dDale et alj|200ll. ICharv & Elbazll200ll and Lagache cold and La- 
gache SB from lLagache et alj2004h and set them as a free parame- 
ter, i.e., for each source, the code chooses the best fitting template 
of all the four libraries. This helps producing more accurate fit- 
ting to the data than choosing only one IR library. The code first 
performs a template-fitting procedure based on a simple \ 2 min- 
imization; then it integrates the resulting best-fit spectra from 8 
up to 1000 /im and directly gives as an output the Lir value. We 
fixed the redshift to the spectroscopic redshift (z s ) of each source 
and made use of multi-wavelength data, including IRAC, Spitzer 
24 /im and recent PEP Herschel data at 100 and 160 (jm. The in- 
clusion of these new data at large wavelengths, which sample the 
critical far-infrared peak of nearby galaxies, allows us to derive the 
Lir with out associated unc e rtainties due to extra polation of the IR 
peak (see lElbaz et"ai1l20ld iNordon et alj|2012l) . We estimate the 
lcr uncertainty of the Lir as half the difference between Lir, sup 
and LiR^nf, where Lm tSU p and Lir^j are the Lir values for 
Ax 2 = 1. In Fig. 1 we show some examples of our SED-fittings. 
Black circles are the observed data, while the best SED fitting are 
the blue and the red lines for the stellar and starburst component, 
respectively. Notice the importance of the PEP data at 100 and 160 
Hm to accurately sample the IR peak of the spectra emission. 



3.2 Stellar Masses 

The galaxy stellar masses have been de rived as explai ned in 
iDominguez Sanchez etld] d201ll) (see also lllbert et alj|2010l) . We 
have used the Lephare code to fit our multi-waveleng th data (from 
NUV to 5.8 /im) with a set of SED templates from lMaraston] d2005l) 
with star formation histories exponentially declining with time as 
SFR oc er t/T . We used 9 different values of r (0.1, 0.3, 1.0, 
2.0, 3.0, 5.0, 10.0, 15.0 and 30.0 Gyr) with 221 steps in age. The 
metallicity assumed is solar. Dust extinction from SED -fitting was 
applied to each galaxy using the ICalzetti et all (2000) extinction 
law. In this paper we allow a maximum E(B — V) value of 1.4. 
This value was chosen as the highest E(B — V) that the code was 
able to handle, given the large grid of parameters used in the SED- 
fitting, and is significantly larger than the highest E(B — V) value 
that we derive for our sample (E(B — V) m ax=0.&5, see Sect. |4]l. 
We imposed to the derived age of the galaxies to be less than the 
age of the Universe at that redshift and greater than 10 8 years. 
The latter requirement avoids having best fit SEDs with unreal- 
istic ally high s pecifi c star formation rates, SSFR — SFR/M. 
The lMarastor] d2005) models include an accurate treatment of the 
thermally pulsing asymptotic giant branch (TP-AGB) phase, which 



1 PEP images and catalogs will be released at 
http://www.mpe.mpg.de/ir/Research/PEP/public_data_releases.php 
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Min 


1 st quartile 


Median 


3 rd quartile 


Max 




0.063 


0.195 


0.267 


0.347 


0.460 


log Mass [Mq] 


8.61 


9.86 


10.10 


10.50 


11.50 


loe SFR (Ha) rMm vr _1 l 


-0.83 


0.42 


0.85 


1.15 


2.32 


log SFR(Lj/{) [Mq yr _1 ] 


-1.03 


0.45 


0.78 


0.99 


1.98 


logSSFRfyr- 1 ] 


-10.80 


-10.10 


-9.65 


-9.36 


-8.31 


Z [log (O/H) + 12] 


8.10 


8.64 


8.76 


8.88 


9.17 


Ha EW [A] 


4.35 


19.60 


31.00 


46.00 


204.00 


H/3 EW [A] 


0.00 


2.35 


4.82 


7.36 


28.30 



Table 1. Properties of the 474 sample galaxies, indicated by their minimum and maximum value, as well as the median, the I s * and the 3 quartile. 
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a width of 1 arcsec, therefore an aperture correction for slit losses 
is applied to each source. Each zCOSMOS spectrum is convolved 
with the ACSI (814) filter and then this magni tude is compared wit h 
the /-band magnitude of the GIM2D fits of ISargent et al.1 j2007l) . 
The difference between the two magnitudes gives the aperture cor- 
rection for each spectrum. This correction assumes that the emis- 
sion line fluxes and the /-band continuum suffer equal slit losses. 

A major problem when deriving SFR(Ha) is the effect of dust 
extinction. Star formation normally takes place in dense and dusty 
molecular clouds regions, so a significant fraction of the emitted 
light from young stars is absorbed by the dust and re-emitted at IR 
wavelengths. A method to measure the extinction is to compare 
the observed ratio of the Ha and H/3 emission lines with the 
theoretical value (R t h = Ha/ H/3 = 2.86). Then the reddening 
toward s the nebular regions E(B — V), following ICalzetti et al.l 
d2000h extinction law, can be written as: 



10.0 100.0 1000.0 
X[^m] 



Figure 1. Example of SED-fitting for 4 sources. The observed SED of each 
source (black filled circles) is shown with the corresponding best-fit solution 
(blue solid line for the stellar part, red solid line for the IR), as well as the 
morphological type (see Sect. 16.61 , the z s and the derived L jjj . 



has a high impact on modeling the templates at ages in the range 
0-3 < t < 2 Gyr, where the fuel consumption in this phase is 
maximum, specially for the near-IR part; although their accurate- 
ness is still debated by a number of rec ent observational studies 
dKriek et alJl201(llMelbourne et al.ll2012h . 



A v = [2.5i? v - x log w (Robs/Rth)]/(kl - k2) 
kl = 2.659 x [-2.156 + (1.509/%) - 
(0.198/A?) + (0.011/A?)] + R v 
k2 = 2.659 x [-1.857 + (1.040/%)] + Rv 
Ai = 0.48613 iim 
A 2 = 0.65628 iim 
Rv = 4.05 
E(B-V) = Av/4.05 



4 SFR FROM Ha EMISSION LINE 

In this paper we will derive the SFR from the Ha emission line 
luminosity corrected for extinction [SFR( Ha) hereafter]. We will 
use the conversion factor derived by iKennicutt et al. I i 1994 and 
lMadauetalJdl998h to obtain our SFR from the Ha emission line: 



SFRiMeyear- 1 ) = 7.9 x 10~ 42 L(H a )(erg s" 1 ) (1) 

Ha emission line fluxes were measured u sing the automatic 
routine Platefit VIMOS {Lamareille etai]|2009l). After removing a 
stellar component using iBruzual & Charlotl d2003"l) models, Plate- 
fit VIMOS performs a simultaneous Gaussian fit of all emission 
lines using a Gaussian profile. The slits in the VIMOS masks have 



Note that this Av and E(B — V) values refer to the attenua- 
tion and reddening towards the nebular regions (as they are derived 
from the nebular emission lines), and not to the attenuation and 
reddening towards the stellar continuum usually found in literature 
and whi ch do not necessari ly take the same values (see Calzetti et 
al. 2000, IWuvtsetai]|201ll) . 

Making use of the spectra of our sources from the 20k survey 
we can measure the observed ratio between the emission lines and 
obtain the extinction values. However, the quality of the spectra is 
not high enough to individually measure the H/3 line for the whole 
sample. To improve the S/N of the spectra we construct average 
spectra by adding spectra of galaxies with similar L/_r values. We 
constructed 6 average spectra composed of 3 1 galaxies each with 
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Figure 2. The black line represents the average spectrum for the lowest 
average hjR before subtracting the stellar component (red), while the spec- 
trum obtained when removing the stellar continuum is plotted in blue. 
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different average L IR values: log L/ H (L )=[9.8, 10.1, 10.4, 10.7, 
11.0, 11.3]. 

For each average spectrum we have measured the Ha and H/3 
fluxes using IRAF, aft er subtracting the stellar c omponent making 
use of more than 230 iBruzual & Chariot! d2003h template models 
with different values of extinction. Fig. [2] shows one of the aver- 
age spectra (black line) before subtracting the stellar continuum 
(red line) and after subtraction of the continuum (blue line). Once 
the stellar component is subtracted from the observed spectra, the 
emission lines can be measured without contamination by the stel- 
lar absorption. Fig. [3] shows the 6 average spectra of our sources 
after performing the continuum subtraction, as well as the values 
of the observed ratio between Ha and H/3 emission lines, Ay and 
E(B — V). 

There is a very clear trend between the extinction and the Lir, 
as can be seen in Fig.[4] where we plot the E(B — V) value derived 
from the average spectra with respect to the average log L/jj of 
each average spectrum. The relation b etween dust extinction an d 
Li r has also been previously studied bv lWang & Heckman| [l996). 
where the authors conclude that the H/3/Ha ratio decreases with 
L/fl (meaning that the dust extinction increases with L/n, in agree- 
ment with Fig. [4j . We fitted the data of Fig. [4] with a straight line, 
obtaining a relation between log Lir and extinction of the form: 



Figure 3. Median spectra in 6 Lir bins. Also shown are the values of R oos , 
Ay and E(B - V). 




log L (L ) 



Figure 4. E(B — V) versus log Ljjj for the 6 average spectra. The thick 
line is the best fit to the data, with a slope m=0.29 (see Eq.[2). 



E(B - V) = 0.29 x logL IR (L 3 ) - 2.54 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 5. Difference between the E(B — V) value derived from the linear 
relation from Eq.|2]and the E(B — V) value measured for each source for 
a control sample of high S/N galaxies in the continuum both in the Ha and 
H/3 spectral ranges. The A(E(B-V)) values are consistent with a gaussian 
distribution with a = 0.20. We assume this dispersion to be the error in 
E(B — V) that we introduce when using an average extinction instead of 
an individual extinction. 

The data follow the derived relation with almost no scatter. We 
can interpolate these data and derive an E(B — V) value for each 
galaxy fr om its Lir and correc t the Ha flux of each galaxy making 
use of the lCalzetti etai] |2000) extinction law: 

F corr (H a ) = F(H a ) x 10 [o.4x3.327 B ( fl -v)] (3) 

We then derive the Ha luminosities using the measured spec- 
troscopic redshift z s of each source, and finally we obtain the 
SFR(Ha) following Eq. 1. 

Our main source of SFR(Ha) error is likely to be the dust 
extinction correction. To test the effect of using a linear relation be- 
tween the extinction and log L/j? instead of using single values of 
extinction derived for each galaxy, we have visually inspected the 
spectra of our sources and we have selected a small sub-sample of 
29 sources with high S/N continuum in both the Ha (S/N c > 2) 
and H/3 (S/N c > 5) spectral ranges, so that we can rely on the mea- 
sured value of the Ha/H/3 ratio to derive E(B — V) for each source 
(Robs,ind)- In Fig.[5]we show the difference between the E(B — V) 
value obtained through the linear relation between log Lir and 
E(B — V) and the E(B — V) value measured for each source 
for this control sample of 29 galaxies. The dispersion is a — 0.20, 
with a small, not significant offset ((A(E(B-V)> = -0.07). This 
small offset confirms that, when using the linear relation between 
extinction and log L/_r instead of using individually measured ex- 
tinctions, we are not introducing a systematic bias but only increase 
the dispersion of the data. We interpret the value a = 0.20 as the 
E(B — V) uncertainty that we introduce in the SFR(Ha) deriva- 
tion. 



5 SFR COMPARISON 

In this section we will compare the SFR(Ha) with the SFR derived 
from the Lir. We have converted our Lir values into SFR mak- 



ing use of the K98 relation, which simply multiplies the L/ij by a 
constant: 

SFRiMoyear' 1 ) = 4.5 x 10" 44 L/ fl (er ff s _1 ) (4) 

In Fig. [6] we compare the SFR(Ha) with the SFR(L/ fl ). The 
blue line is the one to one relation, while the red line is the best 
fit to our data. For the clarity of the figure we do not include the 
errors for each source, but show the median error on log SFR(L/.r) 
and log SFR(Ha) at the right hand of the plot. We have performed 
a linear fit taking into account the errors in both SFR(L/.r) and in 
SFR(Ha). The latter include also the uncertainty introduced when 
using median extinctions (as explained in Sect. [4}- The slope of our 
relation is m — 1.01 ± 0.03. Also the normalizations of the two re- 
lations are well consistent with each other with a very small offset 
of ~ —0.08 ± 0.03 dex. Therefore, there is an excellent agree- 
ment between the two SFR indicators. It is interesting to notice 
the relatively small dispersion of the data (a — 0.28), compara- 
ble with the uncertainties derived for the SFR(Ha) (with a mean 
error ~ 0.30). This is mainly due to the large uncertainties intro- 
duced in the SFR(Ha) derivation when using a linear relation be- 
tween log L/ r and the dust extinction values instead of using the 
E(B — V) values for each galaxy. Considering the wide range of 
SFRs studied (about two orders of magnitude) and the systematic 
errors that could in principle affect our sample, we confirm that the 
two SFR estimates are in an excellent agreement, even if both of 
them are based on very simple recipes that do not take into account 
the galaxies intrinsic properties, but are a mere transformation from 
luminosity to SFR. 

Our results a re consistent with the previous work by 
lKewlevetalJl2002h . where the authors compare SFRs from Ha lu- 
minosity and Ljh for a local sample of galaxies from the Nearby 
Field Galaxy Survey. After reddening correction, they derive a 
slope of 1.07 ±0.03 and a normalization constant of —0.04 ±0.02, 
in a very good agreement with the values that we have obtained. 

It should be mentioned that, since the dust-corrected SFR(Ha) 
is, in the absence of measurement uncertainties, per definition, the 
total SFR, while the SFR(L/_r) is only the obscured part of the 
total, the good one-to-one correspondence between SFR(Ha) and 
SFR(L/i?) implies that the unobscured SFR is much lower than the 
obscured SFR. This is true since our sample is PEP-d etected and 
is in agreement with the result of ICaputi et al.l d2008l) . where the 
authors show that the SFR form the UV and the total SFR (IR+UV) 
differ by a factor ~ 10 on average for a 24/im selected sample from 
the zCOSMOS-bright 10k survey. When studying a random, not IR 
detected sample of galaxies, the emission from young stars coming 
out unobscured in the UV may be significant, in which case the 
SFR(L/_r) would not be in such a good agreement with the total 
SFR(Ha). 

We have considered the possible bias in our results which 
might be induced by a selection of the sample in the far-IR bands. 
For the sources with a detected Ha flux, but not detected in the 
far-IR, we have derived the SFR from the Ha luminosity and then 
transformed this into Lir assuming that the K98 relation is correct 
(i.e. inverting Eq.@. There is a quite robust and well defined corre- 
lation between Lir and 100 and 160 /im fluxes for a given redshift. 
Therefore, knowing the redshift and the Lir, we have calculated 
the approximate PEP fluxes that these sources should have if they 
would follow the K98 relation. 99% of the sources detected in Ha 
but not detected in the far-IR have estimated PEP fluxes lower than 
the PEP detection limit, i.e. they are consistent with following the 
K98 relation and being undetected in the far-IR. 
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Figure 6. SFR(L/j{) vs SFR(Ha).The blue line represents the one to one relation, while the red line is the best-fit to our data. Also shown in the plot the slope 
of the best fit (m = 1.01 i 0.03) and the dispersion of the relation (a = 0.28). The median error on log SFR(L/^j) and log SFR(Hct) is shown towards the 
lower right of the plot. 



We have also studied a possible aperture correction effect. The 
IR fluxes at 100 and 160 /im from Herschel are integrated fluxes 
over the whole galaxy due to the large PSF. Instead, as explained 
in Sect. 4, the Ha is measured with a slit of 1 arcsec width and 
then aperture corrected. This correction assumes that the radial Ha 
flux distribution is the same as the radial distribution of the con- 
tinuum in the 1(814) ACS for each galaxy. However we know that 
galaxies have internal structures, specially affecting the star form- 
ing regions where the Ha emission line is observed. Larger galaxies 
should be more affected by aperture correction as the observed Ha 
flux is measured in a smaller part of the galaxy, while for small 
galaxies almost all the light enters into the slit. We have tested 
the aperture correction effect by comparing the difference between 
SFR(L/_r) and SFR(Ha) with the angular sizes of the galaxies. We 
did not see any systematic effect, i.e. the A(SFR) does not depend 
on the galaxy sizes. We therefore conclude that the aperture cor- 
rection does not systematically affect the slope of the SFR(L/_r)- 
SFR(Ha) relation. 

Finally, we remind that our sample misses about 25% objects 
with PACS detection but Ha S/N < 3. These may include 'aging' 



objects with low SSFR, extending trends discussed in Sect. 6.5, as 
well as some of the most dust obscured systems. 



6 SFR(L /fl )-SFR(Ha) RELATION DEPENDENCES 

We have seen that both SFRs estimates agree very well with each 
other for the sample of galaxies as a whole. In this section we aim 
to detect which intrinsic properties of the galaxies could affect the 
SFR(Ha)-SFR(L/ fl ) relation. 

6.1 Redshift 

In Fig.|7](upper panel) we show log SFR(L/_r) vs log SFR(Ha) for 
different redshifts, where the color code is explained in the plot. To 
better appreciate the effect of z, we plot in the lower panel the ratio 
between log SFR(L/_r) and log SFR(Ha) versus z. The red dots 
are the running mean of the A(SFR)=log SFR(Lr fl )-log SFR(Ha) 
every 20 values of z. The error bar on the left hand of the plot repre- 
sents the typical mean error of the running mean (a rm , z = 0.06). 
There seems to be no significant dependence of the relation on the 
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Figure 7. Upper panel: log SFR(L/ fi ) vs log SFR(Hct) for different redshift 
bins. Dashed and continuous lines represent the best fit to our data and the 
one to one relation, respectively. Lower panel: A(SFR) =log SFR(L//j)- 
log SFR(Ha) versus z. Red dots are the running mean of A(SFR) every 
20 values of z. Also shown at the left hand of the plot is the typical mean 
error of the running mean. 

redshift, but a random scatter around zero. We conclude that the 
SFR(L/_r)-SFR(Hq) relation is not highly affected by the redshift 
of the galaxies, at least up to z = 0.46, which is the maximum 
considered z of our sample. 

6.2 Mass 

In Fig. [8] we show the ratio between log SFR(L/h) and log 
SFR(Ha) versus the stellar mass. The black dots are the running 
mean of the A(SFR)=log SFR(L 7fl )-log SFR(Ha) every 20 values 
of mass and the error bar on the left hand of the plot represents 
the typical mean error of the running mean (oy m , m = 0.06). In 
this case, there is not only a random scatter around zero, but there 
is a slight trend with mass, in the sense that for low mass sources 
the SFR(L/_r) is lower than the SFR(Ha). We have performed a 
Kolmogorov-Smirnov test to the distribution of the ratio between 
the two SFR indicators for the low mass (log M < 10) and high 
mass (log Mass > 10) galaxies, finding that the probability P that 
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Figure 8. log SFR(Ha)-log SFR(Lffj) versus log Mass. The black dots are 
the running mean of A(SFR) every 20 values of stellar mass, while the 
colors represent different metallicity values. Also shown at the left hand of 
the plot is the typical mean error of the running mean. 

the two distributions come from the same population is P — 0.037, 
meaning that the behavior for the low and high mass galaxies is dif- 
ferent at 2.1cr. 

We should remind that our SFR(Ha) is also affected by errors 
and sy stematic uncertainties on its derivation. iBrinchmann et al.l 
d2004l) (B04 hereafter) have shown that, differently from what is 
assumed in Eq.[4] the ratio of observed Ha luminosity to SFR, i.e. 
the efficiency r\ — Lhc/ SFR, is not constant, but depends on 
the mass of the galaxy, with higher values for lower mass systems. 
Therefore, the mass dependence that we observe could be due to 
the fact that we are over-estimating the SFR(Ha) for low mass sys- 
tems. 

To address this issue, we have re-calculated our SFR(Ha) 
making use of the B04 recipes. In this work, the authors derive 
different 77 values for different mass ranges. Even using mass- 
dependent efficiencies, the disagreement between the two SFR in- 
dicators at lower masses is still present, while the slope of the 
SFR(L/_R)-SFR(Ha) relation remains consistent with the one to 
one relation (m = 0.95 ± 0.03). As a further improvement, in 
B04 the authors also derive a mass dependent extinction value for 
the Ha emission line. We have also calculated the SFR values using 
both the mass dependent efficiency and the mass dependent extinc- 
tion from B04 (instead of using the extinction values derived from 
Eq.[2}- When including the extinction values derived from B04, the 
disagreement at lower masses is slightly reduced, but now the val- 
ues of the SFR(L/_r) are significantly higher than the SFR(Ha) at 
higher masses. Moreover, the slope of the SFR(L/_R)-SFR(Ha) re- 
lation derived using the B04 extinction values is significantly larger 
than 1 (m = 1.26 ± 0.05). This is mainly due to the fact that the 
B04 extinction values at high Lir are lower than the extinction 
values that we derive from Eq.|2] 



6.3 Metallicity 

A physical parameter related to the mass is the metallicity. Metal- 
licity an d mass are correlated by the well known mass-metallicity 
relation dTremonti et al.|[2004h . with more massive galaxies being 
more metal-rich, at least up to M~ 10M©. Above this mass, metal- 
licity still increases with mass, but the slope of the relation is sig- 
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Figure 9. log SFR(L/ft)-log (Ha) vs metallicity. Red dots are the running 
mean of A(SFR) every 20 values of Z. The colored lines represent the 
result for the chemo-spectrophotometric model for three spiral galaxies of 
different masses (see Sect. 16. 4| . 



nificantly flatter than at lower masses. Therefore, the difference 
between SFR(Ha) and SFR(L/_r) could also be due to the differ- 
ent metal content of the galaxies. To test this effect we have made 
use of the COS MOS metallicit y catalog by Nair et al. (in prepara- 
tion), using the lDenicolo et alj d2002[) method to estimate metallic- 
ity from [Nil] and Ha emission lines. In Fig.[8]we have plotted our 
sources with different colors for different metallicity values (black 
crosses represent galaxies for which the metallicity measurement 
is not reliable, 48 % of the sample). It can be seen that the lower 
values of the SFR(L/_r) at low masses is mainly driven by a small 
number of very metal-poor galaxies, while most of the most metal- 
rich galaxies lie above the zero point. 

To better appreciate the dependence of the two SFR indica- 
tors from metallicity, in Fig. [9] we show the ratio between the two 
SFR indicators versus the metallicity. The red dots are the run- 
ning mean of A(SFR)=log SFR(L/ fl )-log SFR(Ha) eveiy 20 val- 
ues of Z. Also shown in the left hand of the plot is the typical 
mean error of the running mean {o rm , z = 0.07). It can be ob- 
served that for metal-poor galaxies the SFR(L/_r) values are lower 
than the SFR(Hq), while for metal-rich galaxies the SFR(Lrjj) val- 
ues are larger. The average difference between the two SFR indi- 
cators varies by ~ 0.6 dex when moving from the metal-poor to 
the metal-rich galaxies, meaning that the metallicity has a very im- 
portant effect when deriving the SFR. This is consistent with the 
fact th at more metal-rich galaxies usually hav e higher dust to gas 
ratios ( iBaugh et al.ll2005llSchurer et alj|2009h . meaning that more 
metal-rich galaxies are more efficient in absorbing and re-emitting 
the light emitted by young stars at IR wavelengths (at least until 
the optical depth in the UV goes over r > 1). Therefore, for the 
same SFR values, a metal-rich galaxy emits more in the IR; thus 
the SFR derived from the L/_r is over-estimated with respect to 
the SFR from the Ha luminosity. We remind that the K98 relation 
is calibrated for solar metallicity and we note that for that value 
(12+log(0/H)= 8.69. lAsplund et alj[2009l) the agreement between 
the two SFR indicators is very good. 

It could be argued that the observed trend with metallicity 
could be also affected by using the linear relation between log 
L/n and E(B — V) for all the galaxies. For metal-poor galaxies, 



the dust extinction correction that we apply could be larger than 
the actual one; thus the SFR(Hq) that we derive would be biased 
high, leading to the low values of the ratio between the two SFR 
indicators at low metallicities. However, even when using the B04 
recipes to derive the SFR, which include both a mass-dependent 
efficiency and a mass-dependent dust extinction, the ratio between 
the two SFR estimates was still very different for metal-poor and 
metal-rich galaxies. We eliminated the scatter due to the metallic- 
ity dependence on the SFR comparison, finding a reduction of the 
dispersion of 14% (from a = 0.28 to a = 0.24). We conclude that 
metallicity is a key parameter when deriving the SFR and suggest 
that the K98 relation cannot be applied blindly when considering 
very metal-poor/rich galaxies. 



6.4 Comparison with a model for dust evolution in spiral 
galaxies 

A suite of chemical evolution models for spiral galaxies able to 
reproduce the observed evolution of the mass-metallicity relation 
has been presented in Calura et al. (2009). In that work, three spi- 
ral models, representing galaxies of different masses, have been 
designed in order to reproduce a few basic scaling relations of 
disk galaxies across a stellar mass range spanning from 10 9 Mq to 
10 11 Mq. In such models, it was assumed that the baryonic mass 
of any spiral galaxy is dominated by a thin disk of stars and gas in 
analogy with the Milky Way. The disk is approximated by several 
independent rings, 2 kpc wide, without exchange of matter between 
them. The timescale for disk formation is assumed to increase with 
the galactocentric distance, according to the "inside-out" scenario 
(Matteucci & Francois 1989). The models for spiral disks include 
dust production, mostly from core-collapse supernovae and inter- 
mediate mass stars, restoring significant amounts of dust grains dur- 
ing the asymptotic giant branch phase, as well as dust destruction 
in supernova shocks and dust accretion. The pres criptions for dust 
production are the same as in lCalura et al. I d2008l) . where a chemi- 
cal evolution model for dust was presented, able to account for the 
observed dust budget in the Milky Way disk, in local elliptica ls and 
dwarf galaxies. In a subsequent work. ISchurer et al. (2009) com- 
bined the predicted dust masses and the star formation history of 
the MW galaxy with a spectrophotometric evolution code that in- 
cludes dust reprocessing (GRASIL). 

The relation between log [SFR(L/n)/SFR(Ha)] and metallic- 
ity observed in spiral galaxies is compared with the results from 
chemo-spectrophotometric models of three spiral galaxies of dif- 
ferent masses in Fig. [9] In the model, the SFR(L/_r) quantity is 
calculated according to Eq. [4] where the total IR luminosity Lm 
of the spiral galaxy has been calculated by integrating the SED in 
the range 8 — 1000/im. In Fig. [9] the metallicity (expressed as log 
(O/H) +12) is a direct output of the chemical evolution models of 
Calura et al. (2009). Here, the quantity SFR(Ha) is the intrinsic 
SFR ip, calculated by means of the Schmidt-Kennicutt law 



/ « 
ip oc a 



(5) 



where a is the gas density per unit surface and k ~ 1.5. In agree- 
ment with what is observed in real galaxies, in the three models 
the ratio log [SFR(L/_R)/SFR(Ha)] grows with metallicity. This is 
due to the fact that, as galaxies evolve, their interstellar metallic- 
ity rises and their star formation rate decreases, accompanied by a 
milder decrease of the FIR luminosity. In the model, the smaller 
decrease of the FIR luminosity with respect to the SFR reflects the 
fact that the FIR luminosity is only to a first order approximation 
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Figure 10. log SFR(L /fl )-log SFR(Ha) vs SSFR. Colors represent the 
galaxy stellar ages as derived from the SED-fitting. The black dots are the 
running mean each 20 values of SSFR and the error bar at the left hand of 
the plot is the typical mean error of the running mean. 



proportional to the SFR; in fact, significant contributors to dust pro- 
duction and heating are also intermediate mass stars, whose total 
mass within a galaxy reflects the integrated star formation history. 
Our models predict a larger log [SFR(L/_R)/SFR(Ha)] value at a 
fixed metallicity for lower mass galaxies. It is also worth noting that 
the observed increase of the log [SFR(L/_r)/SFR(Hq)]- metallicity 
relation is stronger than the model results, with a substantial num- 
ber of the observed systems presenting log [SFR(Lr_f?)/SFR(Ha)]> 
0.1, not reproduced by our models. The majority of these systems 
have metallicities 12 + log(0/H) > 8.5, but do not present par- 
ticularly high SFR or stellar mass values. There are two possible 
reasons for this discrepancy between the model results and the ob- 
served galaxies. 

One possibility is that the star formation history of the three 
spiral models does not represent the true star formation history of 
the observed galaxies; a significant number of real galaxies may 
have more complicated star formation histories, for instance char- 
acterized by episodic starbursts, as opposed to our models, char- 
acterized by smooth star formation histories. Another possibility 
is that, for these relatively high metallicity systems, the corrected 
SFR(Ha) values may represent lower limits to the intrinsic SFR 
values. It could be that, the larger the metallicity and the dust con- 
tent of the galaxy, the more likely are the star forming regions to 
become thick to Ha radiation. One consequence may be that par- 
ticularly optically thick regions may be missed, and that the total 
SFR rate determined from Ha emission may be underestimated. 
At present, it is very difficult to assess which one of the two pos- 
sibilities may be the main cause of our discrepancy. A useful test 
in this regard would imply the use of spectro-photometric models 
includ ing the effects of d ust extinction on nebular emission lines 
(e.o. lPanuzzo et al1l2003h . The combined use of such models with 
chemical evolution models including dust evolution will be an in- 
teresting subject for future work and will likely shed more light on 
the main subject of this paper. 



6.5 SSFR 

The SSFR is a very interesting parameter as it does not only ac- 
count for the number of stars being formed, but also takes into ac- 
count the efficiency of a galaxy of a given mass to form stars with 



respect to its total mass content. It gives information about the star 
formation activity of a galaxy, allowing us to divide them into star- 
forming, active or intermediate galaxies. It is reasonable to consider 
that galaxies with different star-formation modes may show differ- 
ent luminosity-SFR relations. In Fig.QJJwe show a similar plot to 
Fig.[8]where now the considered parameter is the SSFR and the col- 
ors represent different ages of the galaxies derived by means of the 
SED-fitting (calculated when deriving the galaxy stellar masses as 
explained in Sec. |3.2t . Also shown is the typical mean error of the 
running mean (<J rm ,ssFR = 0.06). The SFR value that we use in 
this plot is the mean value of SFR(Ha) and SFR(L/_r). This value 
has been chosen to avoid a systematic trend with one of the two 
SFR indicators, which are present in the y axis. As expected, all of 
our galaxies have high values of SSFR and would be classified as 
star-forming galaxi es following different criteria recently found in 
litera ture (see e.g. iDo mmguez Sanch ez~et al.ll20lH . IPozzetti et alj 
[201(J). 

We observe values of the SFR(L/_r) larger than the SFR(Ha) 
by ~ 0.2 dex for galaxies with small values of SSFR (< 
10.0 [yr -1 ]). A possible explanation for these higher values of 
SFR(L/_r) at low SSFRs is the fact that the old stellar population 
may significantly contribute to the L of a galaxy as already men - 
tioned by K98 (see also a recent work by da Cunha et alJrioTlh . 
For the galaxies with low SSFR, the heating by stars older than 
10 7 yr in the diffuse inter stellar medium can be comparable to the 
dust luminosity from the absorption of light from young stars; thus 
the value of the SFR(L/_r) is overestimated as it includes emis- 
sion from the old stellar population and not only from the recently 
formed stars. In fact, as it can be seen in Fig.QJj] the oldest galaxies 
(redder colors) populate the low SSFR region of the plot, while the 
galaxies with the youngest stellar populations (bluer colors) have 
always larger SSFR inferred values. Another possible explanaition 
for the observed trend between the two SFR indicators with SSFR 
could be due to spikes in recent SFR. A galaxy undergoing a cur- 
rent burst (high SSFR) will have a high Ha flux compared with its 
Lir and the log [SFR(L/_R)/SFR(Ha)] will take negative values, 
with the opposite trend happening for low SSFRs. Again, this trend 
is not reduced when using B04 recipes to estimate the SFR(Ha). 
Finally, we would like to remark that the results from Figs. [9] and 
[10] are in agreement with the fun damental metallicity relation ob- 
served bv lMannucci et "alll bOlOh . where the authors find that a de- 
crease in metallicity also correlates with an increase in SSFR. The 
SFR(L/_r) seems to be higher, on average, than the SFR(Ha) both 
for galaxies with high metallicity or low SSFR values. 

We conclude that caution must be taken when deriving the 
SFR from the L/_r for galaxies with low values of SSFR, as the 
emission from the older stellar population may overestimate the de- 
rived SFR. The opposite effect may happen for galaxies with recent 
satr formation bursts, i.e. high SSFR values. 



6.6 Morphological type 

Another important property that could affect the SFR indicators is 
the morphological type, as different morphologies are associated 
with different star formation histories and geometries and therefore 
different SEDs and IR emission. To assess the effect of morphology 
on the SFR derivation we have made use of the morphological cata- 
log by Nair et al. (in pr e parati on), who, following a similar scheme 
as in iNair & Abraham! (120 lOh . visually classified the galaxies be- 
longing to the COSMOS field, dividing them into different types 
(Ell, SO, Sa-Sd, Irr). In FigEJJwe show the usual log SFR(L /fl ) vs 
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Figure 11. log SFR(Lrft) vs log SFR(Ha) for sources classified as Sa, Sb, Sc or Sd/Irr. Red dashed and continuous lines are the best fit to our data and the 
one to one relation, respectively. Also shown the derived values for the slope (m) and its error (<r m ), as well as the scatter of the relation (a) and the number 
of sources for each morphological type. 



log SFR(Ha) plot dividing our sample into different morphological 
types. 

The two SFR indicators are in very good agreement with 
each other for the three spiral sub-samples. For each of them 
the SFR(L/h)-SFR(Hq) relation is consistent, in both normaliza- 
tion and slope, with the one-to-one relation. For the Sd, irregular 
sample (Sd/Irr), the formal best-fit slope becomes much steeper 
(m = 1.62 ± 0.43), but is still consistent with 1 at the 1.5<r level, 
because of the reduced statistics of this sub-sample (we have only 
17 sources classified as Sd or irregular galaxies). The dependence 
of the SFR indicato rs with morphology has also been studied by 
iKewlev et al] 120021) . finding that the ratio between SFR(L/h) and 
SFR(Ha) is almost the same for early (SO-Sab) and late (Sb-Irr) 
type galaxies. If real, one possible explanation for the discrep- 
ancy for the Irr/Sd galaxies is their complex morphology. Irregu- 
lar galaxies are clumpier and less homogeneous than earlier spi- 
ral galaxies, thus making it more difficult and uncertain to mea- 
sure the total Ha flux of the galaxy due to aperture corrections. 
We have also investigated the possible effect of metallicity in the 
Sd/Irr sample, as irr egular local galaxies are found to be metal poor 
dTolstov et aT1l2009l) . However, we did not see a different metallic- 



ity distribution for the Irr/Sd sample than for the whole sample of 
galaxies. Despite the fact that the number of Sd/Irr objects is quite 
small, we observe a slight trend of a steepening of the slope of the 
SFR(L/ _r)-SFR(Hq) relation when moving from early to late spiral 
galaxies. For example, if we analyze together Sa and Sb galaxies 
we obtain a slope of m = 0.93 ± 0.04, i.e. almost a 3cr differ- 
ence with respect to the slope for the Sc sample. This difference 
may be indicative of a weak dependence of the SFR indicators with 
morphology. However, a more significant statistical sample would 
be necessary to assess if the observed steepening of the slope with 
morphological type is a real effect or if it is just due to the few 
sources considered for each morphological type, specially for the 
Sd/Irr sample. 



7 SUMMARY AND CONCLUSIONS 

In this paper we have empirically tested the relation between the 
SFR derived from the L,ir and the SFR derived from the Ha emis- 
sion line. We have studied a sample of 474 galaxies at z=0.06-0.46 
with both Ha and IR detection from 20k zCOSMOS and PEP Her- 
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schel at 100 and 160 /jm. The L/_r has been derived by integrating 
from 8 up to 1000 /im the best-fitting SED to our IR data from 
IRAC, Spitzer and Herschel. We have derived our SFR from the 
Ha extinction corrected emission line. We have constructed 6 me- 
dian spectra with different median hm values and we have found 
a very clear trend between E(B — V) and log Lir. This allows us 
to estimate E(B — V) values for each galaxy (see Eq. [2j even if 
the quality of the spectra does not allow to accurately measure the 
H/3 emission line for each galaxy. Our main conclusions are: 

• We have compared the SFR(L/_r) and SFR(Hq) finding an ex- 
cellent agreement between the two SFR estimates for the bulk of 
the studied galaxies, with the slope of the SFR(L/_r)-SFR(Hq) re- 
lation m — 1.01 ± 0.03 and the normalization constant of a = 
—0.08 ± 0.03. This means that the simple recipes used to convert 
luminosity into SFR by simply scaling relations are consistent with 
each other. The agreement between the two SFR estimates also im- 
plies that the assumptions, for example on star formation histories, 
that were made to derive the SFR recipes (Eq.[T]and|4} should not 
be very different from the properties of the studied sample. The 
comparison of the SFR indicators at low redshift is crucial to test 
the validity of the SFR estimates which are often used in high red- 
shift studies. The main result that we obtain is that both methods are 
in very good agreement for the low redshift sample, allowing us to 
extend their validity at higher redshifts. This allows us to derive 
the SFR of distant galaxies when only the Lir or the Ha informa- 
tion are available, with the assurance that we are not introducing 
important systematic effects. 

• There seems to be no dependence of the SFR indicators with 
redshift, at least up to z ~ 0.46, which is the limit of our sample. 

• The stellar mass seems to have a small influence on the 
SFR(L/_r)-SFR(Hq) relation in the sense that for low mass galax- 
ies the SFR(L/ r) values are lower than the SFR(Ha). However, the 
observed dependence of the SFR indicators with the stellar mass 
seems to be mainly driven by metallicity. 

• The metallicity, in fact, seems to be the parameter that most 
influences the SFR comparison, with average values of log 
[SFR(L /fl )/SFR(Ha)] which differ by ~ 0.6 dex from metal-poor 
to metal-rich galaxies. This is due to the higher efficiency in absorb- 
ing and re-emitting the light from young stars in the IR for more 
metal-rich galaxies. The dispersion of the SFR(L/_r)-SFR(Hq:) re- 
lation is reduced by 14% when taking into account the scatter due 
to the metallicity. We stress that caution must be taken when deriv- 
ing SFR for very metal-poor/rich galaxies. 

• We show how the behavior of the observed SFR(L/fl)/SFR(Ha)- 
metallicity relation finds a natural explanation within the frame of 
a complete theoretical model for dust evolution in spiral galaxies. 

• The SFR(L/_r) values are larger than the SFR(Ha) for galaxies 
with low SSFR values (log SSFR < 10 [yr -1 ]) probably due to the 
non negligible contribution to the IR emission from the old stellar 
population. Stars older than 10 7 yr may significantly contribute to 
the hi r in systems with little star formation activity. The difference 
between the two estimates could also be affected by the effect of 
short and intense bursts of star formation, i.e. galaxies with high 
SSFR values. 

• When separately studying the SFRs indicators as a function of 
morphological types, we find that they are in excellent agreement 
with each other for the three spiral sub-samples. For each of them 
the SFR(L/_r)-SFR(Hq) relation is consistent, in both normaliza- 
tion and slope, with the one-to-one relation. For the Sd, irregular 
sample (Sd/Irr), the formal best-fit slope becomes much steeper 



(m = 1.62 ± 0.43), but is still consistent with 1 at the 1.5a level, 
because of the reduced statistics of this sub-sample. 
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